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Abstract 


Mechanisms  of  heat  generation  and  distribution  in  carbon-fiber-based  composites  subjected 
to  an  alternating  magnetic  field  are  considered.  A  model  that  predicts  the  strength  and 
distribution  of  these  heat  sources  in  the  plane  of  the  cross-ply  laminate  configurations  has  been 
developed  and  verified.  In  this  analysis,  the  fibers  in  a  cross-ply  pair  are  treated  as  a  grid  of 
conductive  loops  in  the  plane.  Each  such  conductive  loop  uses  the  alternating  magnetic  field  to 
produce  a  rotational  electromotive  force  that  induces  electric  fields  in  the  polymeric  regions. 
Induced  electromagnetic  energy  is  converted  into  thermal  energy  through  dielectric  losses  in 
polymeric  regions  between  the  carbon  fibers  in  the  adjacent  orthogonal  plies  that  the  conductive 
loops  comprise.  Each  possible  conductive  loop  is  accounted  for,  and  the  resulting  superposition 
of  potential  differences  at  the  nodes  leads  to  the  in-plane  profile  of  the  electric  field  in  the 
polymeric  regions.  Data  from  AS4  graphite-reinforced  polyetheretherketone  (PEEK)  laminate 
surface  temperature  measurements  using  liquid  crystal  sheets  compare  qualitatively  with  the 
theory. 


u 


Table  of  Contents 


Page 

List  of  Figures .  v 

1.  Introduction .  1 

2.  Formulation  of  Planar  Grid  Model . 3 

3.  Model  Predictions .  6 

4.  Convergence  Analysis .  1 1 

5.  Parametric  Analysis .  13 

5.1  Coil-to-Specimen  Area  Ratio .  13 

5.2  Location  of  Coil .  15 

6.  Experimental  Support .  18 

7.  Summary .  21 

8.  References . . .  23 

Distribution  List .  25 

Report  Documentation  Page .  35 


in 


Intentionally  left  blank. 


IV 


List  of  Figures 


Figure  Page 


1.  Induced  Current  Due  to  a  Transverse  (Normal  to  the  Plane)  Magnetic  Field 

in  a  [0/90]s  Cross-Ply .  1 

2.  An  Electrical  Network  Analog  to  a  [0/90]  Cross-Ply  With  a  5  x  5 

Grid  Size  Representation  in  the  Plane .  2 

3.  Comparison  of  Laminate  Configuration  to  a  Representative  Planar  Grid 

in  the  Ply-Ply  Interaction  Submodel .  4 

4.  Process  Flow  Diagram  for  Planar  Grid  Model .  7 

5.  A  5  x  5  Grid  Representation  Showing  the  Distribution  of  Unit  Cell 

Magnetic  Flux .  8 

6.  Columnar  Plot  of  Equation  (2)  for  the  Output  of  the  Planar  Grid  Model’s 

Computer  Code .  9 

7.  Plot  of  Main  Diagonal  Elements  of  Equation  (2)  for  the  Nondimensional  Output 

of  the  Planar  Grid  Model’s  Computer  Code  for  a  5  x  5  Grid  Size .  9 

8.  Prediction  of  the  Points  Highest  Heating  for  a  5  x  5  Grid  Size  Representation 

(Darkened  Circles) .  10 

9.  Superposition  of  Diagonal  Nondimensional  Voltage  Distributions 

for  Various  Grid  Densities  From  Planar  Grid  Model .  12 

10.  Plot  of  Average  Nondimensional  Voltages  for  Various  Grid  Densities  From  the 

Planar  Grid  Model .  12 

1 1 .  Plot  of  Percent  Error  of  Average  Nondimensional  Voltage  for  Various  Grid  Sizes 

From  the  Cocon vergent  Solution  at  Infinite  Grid  Fineness .  13 

12.  Results  of  Study  to  Determine  Minimum  Applicable  Grid  Size  for  the  Situation 
in  Which  the  Coil  Completely  Covers  the  Specimen 

(Coil-to-Specimen  Area  Ratio  of  Unity) .  14 

13.  Predictions  of  the  Effects  of  Changing  the  Size  of  the  Helmholtz-Type  Coil 
With  Respect  to  the  Size  of  the  Laminate  Specimen 

for  Centered-Coil  Experiments . 14 


v 


Figure  Page 

1 4.  Averaged  Results  of  Figure  1 3  Showing  Decrease  in  Average  Nondimensional 

Voltage  With  Increasing  Size  of  Coil  With  Respect  to  the  Specimen .  15 

15.  Predictions  of  Voltage  per  Unit  Magnetic  Flux  for  Varying  the  Location 

of  the  Hemholtz-Type  Coil  on  a  Cross-Ply  Specimen .  16 

1 6.  The  “Accuracy  Zone”  for  Coil-to-Specimen  Size  Ratio  as  it  Relates  to  the 

Grid  Density  Used  in  the  Planar  Grid  Model .  17 

17.  A  10  x  10  APC-2  Tape  Grid . . .  18 

18.  The  Predicted  Nondimensional  Voltage  Profile  From  the  Planar  Grid  Model 
for  the  10  x  10  Grid  Representation  Used  to  Model  the  Tape  Segments 

of  Figure  17 .  *9 

1 9.  Results  of  Liquid-Crystal  Thermal  Measurement  Observations 

for  a  10-cm  Helmholtz  Coil  on  a  20-cm-Square  [0/90]  Cross-Ply  Laminate .  20 


vi 


1.  Introduction 


With  the  advent  of  advanced  thermoplastic-based  composites,  much  research  has  been 
directed  to  take  advantage  of  their  unique  postprocessing  attributes.  Thermoplastic  resins  are 
stable  high  molecular  weight  polymers  that  retain  their  chemical  identity  during  processing. 
Since  the  fully  polymerized  thermoplastic  resins  do  not  form  cross-linked  networks,  they  can  be 
softened  and  reformed.  Alternating  magnetic  fields  provide  a  localized,  noncontact,  and 
expedient  source  of  heating. 

It  has  been  established  [1]  that  the  primary  mechanism  of  heating  in  continuous-fiber 
laminated  systems,  such  as  AS4/polyetheretherketone  (PEEK)  carbon/thermoplastic,  is  primarily 
due  to  dielectric  losses  in  the  polymer.  This  was  shown  to  be  true  as  long  as  dielectric 
breakdown  did  not  occur  in  the  polymer.  This  “local  theory”  of  heating  established  the 
mechanism  by  which  electromagnetic  energy  is  converted  into  heating  in  the  locality  of  the 
fiber-fiber  intersection,  as  shown  in  Figure  1. 


Figure  1.  Induced  Current  Due  to  a  Transverse  (Normal  to  the  Plane)  Magnetic  Field  in  a 
[0/90]s  Cross-Ply.  A  Matrix-Rich  Region  of  Thickness  h  Exists  Between  the 
orthogonal  plies  in  such  a  laminate. 
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The  local  theory  provides  the  basis  for  a  “global  model”  of  heat  generation  in 
continuous-fiber  cross-ply  laminated  composite  systems.  A  global  model  is  needed  to  determine 
the  value  of  the  alternating  electric  field  across  the  fiber-fiber  intersection  in  all  interfiber 
polymeric  regions  in  order  to  quantify  the  distribution  of  overall  heating  in  the  specimen.  The 
global  model  [2]  consists  of  three  additional  independent  submodels:  (1)  a  fiber  layer  submodel 
to  analyze  the  through-thickness  electromagnetic  response  in  the  composite,  (2)  a  thermal 
submodel  to  determine  the  surface  transient  thermal  response,  and  (3)  the  planar  grid  submodel 
presented  here  to  describe  the  in-plane  response. 

In  order  to  correlate  the  local  theory  of  heating  occurring  due  to  electrical  potential 
differences  between  intersecting  fibers  with  a  compatible  model  of  thermal  generation  in  the 
plane,  the  laminate  is  modeled  as  an  electrical  network  of  intersecting  conductors  with  some 
reasonable  mesh  size.  For  example,  Figure  2  shows  the  electrical  network  analog  associated 
with  a  5  x  5  grid  size  representation.  The  objective  of  this  work  is  to  characterize  the  interaction 
between  adjacent  orthogonal  or  off-axis  plies  in  the  composite  laminate  subjected  to  a  transverse 
alternating  magnetic  field.  The  model  developed  to  determine  this  two-dimensional  ply-ply 
interaction  is  termed  the  planar  grid  model  to  describe  its  use  of  a  finite  grid  to  represent  the 
plane  of  the  laminate  specimen.  The  interaction  between  individual  fibers  through  the  thickness 
of  the  laminate  is  reported  elsewhere. 
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Figure  2.  An  Electrical  Network  Analog  to  a  [0/90]  Cross-Ply  With  a  5  x  5  Grid  Size 

Representation  in  the  Plane.  Note  That  Fiber-Fiber  Intersection  Resistances  Are 
Considered  to  Be  Equivalent  in  the  Ply-Ply  Interaction  Submodel. 
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2.  Formulation  of  Planar  Grid  Model 


Several  effects  are  considered  in  this  analysis.  These  include 

(1)  cancellation  of  electric  fields  in  internal  loops, 

(2)  determination  of  least-resistive  path  with  respect  to  matrix-rich  intersections, 

(3)  determination  of  least-resistive  path  with  respect  to  fibers,  and 

(4)  incorporation  of  current  density  effects 

These  items,  superimposed,  provide  a  view  of  the  planar  heating  pattern  through  the  distribution 
of  the  electric  fields  along  the  various  conductive  paths.  The  planar  grid  model  incorporates 
items  1,  2,  and  3.  Item  4  would  account  for  the  possibility  of  parallel  fibers  within  the  same  ply 
randomly  coming  in  contact  so  that  current  would  have  the  option  of  taking  several  paths  in 
accordance  with  the  effective  resistances  of  the  various  paths.  Such  effects  would  only  perturb 
the  distribution  of  electric  fields  within  a  few  fiber  diameters.  Accordingly,  these  effects  are  not 
considered.  As  a  consequence  of  this  simplification,  a  relatively  coarse  grid  can  provide 
reasonable  estimates. 

The  cancellation  of  linear  electric  fields  in  internal  conductive  loops  is  a  key  element  of  the 
model.  The  applied  alternating  magnetic  field  induces  a  rotational  electric  potential  field  about 
each  possible  conductive  loop,  regardless  of  the  loop’s  effective  resistance.  A  square  grid,  such 
as  that  in  Figure  2  and  in  cross-ply  laminates  illustrated  in  Figure  3,  can  be  divided  into  many 
possible  conductive  loops  of  various  shapes  and  sizes.  The  minimum  number  of  intersections, 
however,  is  four.  Three-sided  paths  are  not  possible  since  two  interacting  unidirectional  plies 
can  only  form  four-sided  and  greater  paths  when  viewed  normal  to  the  plane.  In  consideration  of 
item  2,  the  least-resistive  path  will  generally  be  a  path  consisting  of  the  minimum  number  of 
intersections.  The  resistance  of  the  fiber  lengths  has  been  shown  to  be  negligible  compared  to 
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Figure  3.  Comparison  of  Laminate  Configuration  to  a  Representative  Planar  Grid  in  the 
Ply-Ply  Interaction  Submodel. 


the  resistance  of  the  intersections  so  that  the  lengths  of  the  paths  traveled  are  inconsequential 
(item  3)  when  considering  the  path  resistance. 


An  algorithm  has  been  developed  that  incorporates  all  possible  four-cornered  loops  in  any 
given  grid  and  calculates  and  superimposes  the  induced  potentials.  Figure  3  shows  the  actual 
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laminate  configuration  considered  and  a  representative  planar  grid  model.  Note  that,  in  the 
laminate,  there  are  approximately  24,000  fiber  rows  in  each  20.3-cm-wide  (8  in)  ply,  which 
combine  to  make  up  approximately  8  x  1016  loops.  The  model  assumes  an  n  x  n  grid  (where  n  is 
some  small  integer)  consisting  of  1/4  n2(n  +  l)2  possible  four-cornered  loops.  The  model  of 
Figure  3  illustrates  a  3  x  3  grid  representation,  containing  36  possible  four-cornered  loops.  For 
each  possible  loop  (which  consists  of  four  fiber  lengths  and  four  matrix  regions),  the  following 
operations  are  performed. 


(1)  Calculate  the  planar  rotational  electric  potential  field  (emf)  from  Faraday’s  law: 


where  is  the  area  of  the  loop  and  the  time  derivative  of  magnetic  field  vector  B  is  the 
product  of  the  angular  frequency  and  the  scalar  B. 

(2)  Convert  the  rotational  emf  for  each  loop  to  a  directional  electric  field  vector  along  each 
fiber  length,  which  comprises  the  loop,  by  dividing  the  emf  by  the  loop  perimeter  in 
accordance  with 


e  =  fE-dI. 

(3)  Sum  the  electric  fields  from  all  loops  for  each  fiber  segment  (grid  element)  obtained 
from  each  loop  calculation  (steps  1  and  2). 

(4)  Calculate  the  alternating  potential  differences  across  each  node  in  the  plane. 

Step  4  provides  the  “nodal”  potential  difference  between  fibers  in  adjacent  plies.  A  separate 
model  is  needed  to  determine  how  the  “layers”  of  fibers  through  the  thickness  interact  with  their 
counterparts  in  the  adjacent  orthogonal  (or  off  axis)  ply.  This  through-thickness  model  is 
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described  in  Fink  [2].  The  planar  grid  model  can,  however,  be  further  analyzed  with  the 
realization  that  it  predicts  the  qualitative  nature  of  heating  in  the  plane  of  the  laminate  since 
heating  through  dielectric  losses  [1]  is  directly  proportional  to  the  square  of  the  potential 
difference: 


6  V2 

Wj  a) 

h 

where  Wj  is  the  heat  generation  at  some  node  j;  f$j  is  a  function  of  several  material, 
environmental,  and  microstructural  properties  at  node  j;  V.  is  the  potential  difference  between  the 
fibers  at  node  j;  and  h  is  the  distance  through  which  the  electric  field  created  by  V  exists,  as 
defined  in  Figure  1. 

Although  the  voltages  cannot  be  directly  measured  (due  to  the  high  frequency)  or  the 
existence  of  the  electric  fields  directly  proven,  their  manifestation  as  surface  temperature 
gradients  can  be  observed.  Parametric  studies  were  performed  verifying  the  convergence  of  the 
grid  size  to  low  n  values  at  various  coil-to-specimen  size  ratios  and  coil  locations.  Experimental 
studies  verify  the  location  of  thermal  extremum  in  the  plane,  as  predicted  by  the  algorithm. 

3.  Model  Predictions 

Figure  4  shows  an  outline  of  a  computer  program,  which  performs  the  operations  described 
previously.  Data  representing  the  input  magnetic  flux  from  the  coil  through  each  smallest  unit 
loop  in  the  grid  (each  element)  are  input  to  the  algorithm.  Equation  (2)  is  an  example  input 
matrix  representing  the  5  x  5  grid  of  Figure  2,  with  a  centered  coil  superimposed.  Since  the 
Helmholtz-type  coil  that  was  used  in  the  experimental  work  provides  a  uniform  distribution  of 
flux  in  the  plane,  the  contribution  of  total  flux  to  each  grid  loop  or  element  can  be  calculated. 
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Specimen  and  Coil  Properties; 

*  Ply  Sequence 

*  Size  of  Specimen 

*  Shape  of  Coil 

*  Relation  of  Coil  to  Specimen 
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Magnetic  Reid 
Strength  into  ^ 
Grid  Elements 


Input  Magnetic  Field 
Strength  Distribution  of 
Coil 


Output 
2-D  Potential 
Difference 
Distribution  for  Ply- 
Ply  Interaction 


Convert  Electric 
field  vectors  to 
potential 
differences  at 
nodes. 
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Determine  all 
possible  conductive 
loops  in  grid 
network,  calculate 
electric  field 
vectors,  and  sum  at 
element  edges. 


Figure  4.  Process  Flow  Diagram  for  Planar  Grid  Model. 


Figure  5  shows  the  5  x  5  grid  representation  superimposed  on  a  square  laminated  cross-ply 
specimen  with  a  Helmholtz-type  coil  centered  on  the  specimen.  The  placement  of  the  coil 
determines  the  area  within  which  the  alternating  magnetic  flux  acts  normal  to  the  surface.  If 
each  element  of  the  modelled  grid  has  a  unit  area,  those  elements  that  are  completely  enclosed  by 
the  magnetic  flux  field  [e.g.,  element  (3,3)  in  Figure  5]  are  considered  to  have  a  unit  flux.  Other 
cells  may  be  only  partially  influenced  by  the  flux  field  and  have  proportionate  fractions  of  the 
unit  flux  assigned  to  them  in  amounts  equivalent  to  the  proportionate  fraction  of  element  area 
covered  by  the  homogeneous  flux  field.  Figure  5  shows  these  fractions  for  the  example  under 
consideration. 


Next,  the  algorithm  normalizes  these  values  so  that  the  total  imposed  magnetic  flux  is  a  unit 
value.  These  values  are  then  used  as  input  as  shown  in  equation  (2): 
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0.036 
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0.036 

0 
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0 

0 

0 
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(2) 
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Figure  5.  A  5  x  5  Grid  Representation  Showing  the  Distribution  of  Unit  Cell  Magnetic 
Flux.  Note  That  Element  (3,3)  Has  a  total  Flux  Input  of  Unity. 

Note  that  the  sum  of  all  the  elements  in  the  input  matrix,  equation  (2),  is  unity.  The  code  returns 
the  output  (nodal  voltage  per  unit  magnetic  flux)  per  equation  (3)  and  is  displayed  in  Figures  6 
and  7. 

Note  in  equation  (2)  that  the  coil  was  centered  on  the  specimen  so  that  the  diagonal  plot  of 
Figure  7  provides  an  avenue  for  comparison  with  other  centered-coil  (i.e.,  symmetrical  coil 
placement)  examples: 
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Figure  6.  Columnar  Plot  of  Equation  (2)  for  the  Output  of  the  Planar  Grid  Model’s 
Computer  Code.  Each  Column  Represents  the  Voltage  Between  Plies  for  a 
Planar  Grid  Node  in  the  Plane  of  the  Laminate.  The  Relative  Position  of  the 
Coil  Is  Shown. 


Voltage 
per  unit 
Magnetic  Flux 


Position  Along  Diagonal 

Figure  7.  Plot  of  Main  Diagonal  Elements  of  Equation  (2)  for  the  Nondimensional  Output 
of  the  Planar  Grid  Model’s  Computer  Code  for  a  5  x  5  Grid  Size. 
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Output  =  [Al,s 


'6.72  9.76  7.49  7.49  9.76  6.72 

9.76  14.18  11.00  11.00  14.18  9.76 

7.49  11.00  6.28  6.28  11.00  7.49 

7.49  11.00  6.28  6.28  11.00  7.49 

9.76  14.18  11.00  11.00  14.18  9.76 

6.72  9.76  7.49  7.49  9.76  6.72 


(3) 


where  vre  is  the  nodal  voltage  in  volts  at  node  (r,s),  (0  is  the  angular  frequency,  and  <{>b  is  the 
magnetic  flux  in  webers.  Each  number  in  equation  (3)  (output)  represents  the  potential 
difference  at  each  node  of  equation  (2)  (input).  Note  that  some  amount  of  voltage  exists  at  each 
node  and  that  the  highest  voltages  occur  at  the  comers  of  the  polygon  formed  by  the  orthogonal 
tangents  to  the  coil  or  flux  region  as  shown  in  Figure  8.  This  distribution  indicates  that  the 
nature  of  the  thermal  response  in  the  laminate  is  dependent  upon  the  size  and  shape  of  the  coil 
and  that  the  model’s  prediction  is  a  function  of  the  grid  dimension  used. 
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Figure  8.  Prediction  of  the  Points  of  Highest  Heating  for  a  5  x  5  Grid  Size  Representation 
(Darkened  Circles).  In  a  Test  Specimen,  the  Points  of  Highest  Heating  Fall  at 
the  Points  of  Intersection  of  the  Tangents  (Dashed  Lines)  to  the  Flux  Region 
(Bold  Circle)  Due  to  the  Nearly  Infinitely  Fine  Grid. 
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4.  Convergence  Analysis 


If  nxxny  is  the  size  of  the  grid  (i.e.,  nx  =  5,  ny  =  5  in  Figure  5),  then  the  total  number  of  loops 
possible  is  1/4  (nx)(ny)(nx  +  l)(ny  +1).  For  a  5  x  5  grid,  this  represents  225  loops,  and,  for  a 
20  x  20  grid,  44,100  loops  must  be  considered.  Therefore,  practical  programing  limitations  on 
grid  size  exist.  For  example,  a  20.3-cm-square  (8  in)  cross-ply  specimen,  such  as  that  used  in 
many  of  our  experiments,  would  require  approximately  a  25,000  x  25,000  grid  for  exact 
representation.  A  convergence  study  was  conducted  to  determine  a  minimum  grid  size  required 
to  achieve  sufficiently  accurate  results.  These  grids  were  then  used  in  further  studies. 

Square  grid  sizes  ranging  from  nx  =  ny  =  3  to  16  were  studied.  For  this  study,  a 
coil-to-specimen  area  ratio  of  14.2%  was  used.  (With  the  coil  centered  on  the  specimen,  the 
percentage  of  area  covered  by  the  coil  and,  thus,  by  the  flux  field  was  14.2%.)  Inputting  a 
standardized  unit  flux,  the  amount  of  flux  through  each  element  could  be  determined  as 
described  earlier.  This  provided  the  input  for  each  case  studied.  Figure  9  shows  the 
superposition  of  diagonal  potential  distributions  for  five  of  the  cases  studied.  The  results  are 
plotted  as  straight  lines  connecting  the  data  points  for  ease  of  reading.  A  quantitative  measure  of 
the  convergence  is  possible  by  comparing  the  volume  under  the  surface  plots  for  each  case.  This 
is  equivalent  to  the  average  of  all  voltage  values  in  the  respective  output  matrices.  Figure  10 
shows  these  values  plotted  against  the  square  grid  size.  Convergence  occurs  rather  quickly,  as 
shown  again  in  Figure  11,  where  the  percent  error  from  the  convergence  value  (8.6  in  this 
example)  is  plotted  against  increasing  grid  size.  Although  this  shows  that  a  square  grid  size  of  8 
could  be  used  with  less  than  a  5%  error,  recall  that  this  result  is  valid  only  when  the  area  of  the 
specimen  surface  covered  by  the  coil  is  equal  to  or  greater  than  14.2%.  Other  considerations, 
such  as  the  minimum  number  of  grid  elements  or  nodes  covered  by  the  coil,  must  be  considered. 
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Figure  9.  Superposition  of  Diagonal  Nondimensional  Voltage  Distributions  for  Various 
Grid  Densities  From  Planar  Grid  Model. 
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Figure  10.  Plot  of  Average  Nondimensional  Voltage  Values  for  Various  Grid  Densities 
From  the  Planar  Grid  Model.  Note  That  the  Horizontal  Axis  Is  Placed  at  the 
Point  of  Convergence  (8.6)  on  the  Vertical  Axis. 
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Figure  11.  Plot  of  Percent  Error  of  Average  Nondimensional  Voltage  for  Various  Grid 
Sizes  From  the  Convergent  Solution  at  Infinite  Grid  Fineness.  Note  That  an 
8x8  Grid  Size  Is  Within  a  5%  Error  Limit  From  the  Convergent  Value 
Determined  From  Figure  10. 

5.  Parametric  Analysis 

5.1  Coil-to-Specimen  Area  Ratio.  At  one  extreme  where  Acou/Aspecimen  is  unity,  the  coil 
completely  covers  the  cross-ply  specimen.  As  usual,  it  is  assumed  that  a  homogenous  magnetic 
field  was  produced  by  the  Helmholtz  coil.  Grid  sizes  of  3  through  9  were  studied  for  this 
situation.  The  percent  differences  (error  values)  are  shown  in  Figure  12.  Note  that  a  grid  size  of 
7x7  falls  within  our  5%  error  standard.  (The  maximum  voltages  for  each  grid  size  remained 
constant.)  The  convergence  of  the  model  at  a  fairly  coarse  grid  size  not  only  makes  calculations 
faster  but  also  validates  the  use  of  the  model  for  representing  the  actual  case  of  fibers  forming  a 
much  finer  mesh. 

Figure  13  shows  the  results  of  applying  a  unit  magnetic  flux  to  various  coil-to-specimen  size 
ratios.  In  each  case,  the  same  total  amount  of  flux  is  input  to  the  specimen  but  the  total 
nondimensional  voltage  is  not  constant.  Figure  14  is  a  plot  of  the  coil-to-specimen  size  ratio 
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Figure  12.  Results  of  Study  to  Determine  Minimum  Applicable  Grid  Size  for  the  Situation 
in  Which  the  Coil  Completely  Covers  the  Specimen  (CoU-to-Specimen  Area 
Ratio  of  Unity). 
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Figure  13.  Predictions  of  the  Effects  of  Changing  the  Size  of  the  Helmholtz-Type  Coil  With 
Respect  to  the  Size  of  the  Laminate  Specimen  for  Centered-Coil  Experiments. 
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Figure  14.  Averaged  Results  of  Figure  13  Showing  Decrease  in  Average  Nondimensional 
Voltage  With  Increasing  Size  of  Coil  With  Respect  to  the  Specimen. 

against  the  average  nondimensional  emf  induced  in  the  specimen.  This  result  shows  that 
increasing  the  area  covered  by  the  coil  without  increasing  the  total  input  flux  decreases  the  total 
resulting  emf  energy.  Conversely,  localization  of  the  flux  in  the  plane  increases  both  the  total 
energy  dissipated  by  the  specimen  and  the  gradients  of  heating  in  the  plane. 

5.2  Location  of  Coil.  Three  parametric  studies  were  performed  for  a  6  x  6  grid  with  a  coil 
that  covered  14.2%  of  the  grid  surface.  The  coil  is  placed  centered,  shifted  to  an  edge,  and 
shifted  to  a  comer,  respectively,  in  the  three  cases.  Figure  15  displays  the  three-dimensional 
(3-D)  columnar  plots  for  the  three  cases  with  their  respective  coil  placements.  Note  that  the 
location  of  the  coil  divides  the  total  grid  into  four  quadrants.  If  the  coil  is  symmetrical,  each 
quadrant  has  the  same  amount  of  energy  induced,  regardless  of  where  the  coil  is  placed.  For 
example,  moving  the  coil  to  a  comer  requires  that  one-fourth  of  the  energy  be  dissipated  in  that 
comer. 

Keeping  the  size  of  the  coil  and  its  energy  constant,  but  moving  it  about  in  the  plane  of  the 
specimen,  changes  both  the  average  induced  voltage  and  the  maximum  voltage  value.  The 
maximum  energy  is  felt  by  the  specimen  when  the  coil  is  centered  in  the  plane,  and  the  minimum 
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Coil  Shifted  to  Corner 


Figure  15.  Predictions  of  Voltage  per  Unit  Magnetic  Flux  for  Varying  the  Location  of  the 
Helmholtz-Type  Coil  on  a  Cross-Ply  Specimen.  The  Relative  Size  and 
Placement  of  the  Coil  Is  Shown  for  Each  Case.  Note  the  Division  of  the  Input 
Flux  Into  Four  Equal  Quadrants. 
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total  energy  is  experienced  when  the  coil  is  moved  to  a  comer.  For  the  slight  shifts  in  coil 
position  shown  in  Figure  15,  the  total  energy  disipated  decreases  4%  and  7.5%  for  the  edge  shift 
and  comer  shift,  respectively.  For  a  situation  in  which  all  the  flux  is  forced  into  the  comer 
element  of  the  6  x  6  grid,  the  decrease  in  energy  dissipated  is  44%.  Note  that  the  maximum 
voltage  is  still  increased  as  the  coil  moves  toward  an  edge  (+16%)  or  comer  (+3%).  For  the 
comer-point-flux  case,  the  increase  is  310%;  however,  this  situation  also  involves  decreasing  the 
size  of  the  coil  with  respect  to  the  specimen.  These  observations  explain  the  “edge  and  comer 
effects”  described  in  the  literature  [3,  4], 

The  convergence  of  the  model  for  square  cross-ply  specimens  was  examined  at  two 
coil-to-specimen  area  ratios:  14.2%  and  100%.  The  results  of  these  studies  (Figures  11  and  12 
respectively)  indicate  minimum  grid  sizes  of  8  x  8  and  7x7,  respectively,  for  errors  of  5%  from 
infinite  grid  fineness.  It  appears  to  be  a  reasonable  assumption  that,  for  any  value  of 
coil-to-specimen  area  ratio  between  14%  and  100%,  the  minimum  square  grid  size  would  fall 
between  7  and  8,  as  shown  in  Figure  16. 


Mnimum  Nominal 
Square  Grid  Size 
for  5%  Error 


Coil-to-Specimen  Area  Ratio 


Figure  16.  The  “Accuracy  Zone”  for  Coil-to-Specimen  Size  Ratio  as  It  Relates  to  the  Grid 
Density  Used  in  the  Planar  Grid  Model.  Note  That  a  Minimum  Grid  Size  of 
8  x  8  Is  Accurate  for  All  Coil-to-Specimen  Area  Ratios. 
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For  coil  sizes  smaller  than  14%,  a  steep  increase  in  minimum  grid  size  is  expected.  For 
example,  consider  a  20.3-cm-square  (8  in)  specimen  with  a  1.3-cm-diameter  (0.5  in)  coil.  The 
coil-to-specimen  area  ratio  is  approximately  0.003.  For  the  coil  placed  at  the  center  of  the 
specimen,  a  16  x  16  grid  size  is  necessary  before  any  changes  in  the  result  occur  since,  for  grid 
sizes  coarser  than  16  x  16,  the  coil  diameter  is  less  than  the  smallest  conductive  loop  in  the 
model.  For  square  cross-ply  specimens  with  centered  coils,  grid  sizes  of  7  or  8  are  sufficiently 
accurate  for  coil-to-specimen  area  ratios  greater  than  14%. 

6.  Experimental  Support 

A  10  x  10  APC-2  tape  grid  was  laid  out  between  plates  of  glass,  and  a  magnetic  field  was 
applied  using  the  Helmholtz  coil;  each  tape  length  was  treated  as  a  conductive  element  in  the 
model.  Figure  17  shows  the  tape  layout,  coil  placement,  and  liquid-crystal  thermal  profile.  The 
small  circles  represent  the  points  of  heating,  as  indicated  by  the  liquid-crystal  sheet  in  the 
40-45°C  range.  The  intensity  of  heating  is  thus  indicated  by  the  size  of  the  dots.  Note  the  four 
points  of  highest  heating  and  the  eight  points  of  second-highest  heating. 

Figure  18  shows  the  3-D  mapping  of  the  model’s  prediction,  which  coincides  qualitatively 
with  the  experimental  observation  of  Figure  17.  Only  the  first  two  sets  of  “highest  heating”  are 
shown.  All  lower  nondimensional  voltages,  representing  values  less  than  60%  of  the  maximum, 
are  omitted  for  clarity. 

A  20.3-cm-square  (8  in)  cross-ply  AS4  graphite/PEEK  [0/90]s  laminated  plate  was  examined 
using  a  10.2-cm-diameter  (4  in)  Helmholtz  coil  placed  at  the  center,  edge,  and  comer  of  the 
specimen.  Figure  19  displays  the  results  of  viewing  40-45°C  liquid-crystal  sheets  during  the 
heating.  The  prediction  of  heating  for  each  case  is  given  in  Figure  15.  A  comparison  of  Figure 
19,  with  predictions  of  Figure  15,  indicates  a  close  correlation  between  the  planar  voltage 
distribution  and  heating  in  the  plane. 
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Figure  17.  A 10  x  10  APC-2  Tape  Grid.  The  Large  Circle  Represents  the  Placement  of 
the  Helmholtz  Coil.  The  Small  Circles  Represent  the  Points  of  Heating  as 
Indicated  by  Liquid  Crystal  Sheet  (40-45°C  Range).  The  Intensity  of  Heating 
Is  Thus  Indicated  by  the  Size  of  the  Dots.  Note  the  Four  Points  of  Highest 
Heating  and  the  Eight  Points  of  Second-Highest  Heating. 


Voltage 
per  unit 
Magnetic  Flux 


Figure  18.  The  Predicted  Nondimensional  Voltage  Profile  From  the  Planar  Grid  Model 
for  the  10  x  10  Grid  Representation  Used  to  Model  the  Tape  Segments  of 
Figure  17.  Only  the  First  Two  Sets  of  Highest  Heating  (Those  in  Excess  of  13) 
Are  Shown.  All  Lower  Voltages  Are  Omitted  Here  for  Clarity.  The  Ring 
Above  the  Graph  Indicates  the  Placement  of  the  Helmholtz  Coil  on  the 
Specimen. 
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Figure  19.  Results  of  Liquid-Crystal  Thermal  Measurement  Observations  for  a  10-cm 
Helmholtz  Coil  on  a  20-cm-square  [0/90]  Cross-Ply  Laminate.  The  Ring 
Indicates  the  Placement  of  the  Coil.  The  Contours  Represent  the  Progression 
(From  Inside  To  Outside)  of  Heating  as  Witnessed  From  the  Liquid  Crystal  in 
the  40-45°C  Range.  Note  the  Points  of  Highest  Heating,  as  Predicted  in 
Figure  15. 
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7.  Summary 


A  model  to  predict  the  distinct  in-plane  response  of  a  continuous-carbon-fiber  thermoplastic 
matrix  cross-ply  laminated  composite  plate  to  an  alternating  magnetic  field  has  been  developed. 
This  model  describes  how  the  transversely  applied  magnetic  field  creates  a  rotational  electrical 
potential  field  that  induces  a  distribution  of  linear  electric  fields  and  nodal  linear  potential 
differences  between  crossing  fibers  in  the  plane  of  the  cross-ply  laminate.  The  planar  grid  model 
is  represented  by  an  algorithm  that  considers  all  the  possible  conductive  loops  in  the  planar 
system  of  crossing  fibers  by  assuming  a  coarse  grid  density.  The  solution  of  this  algorithm  was 
shown  to  converge  at  a  finite  grid  density,  depending  upon  the  size  of  the  coil  with  respect  to  the 
specimen  and  upon  the  spatial  placement  of  the  coil  on  the  specimen.  The  size  and  placement  of 
the  coil  were  also  shown  to  significantly  (and  predictably)  affect  the  strength  and  distribution  of 
the  electromagnetic  response  in  the  plane.  This  response  was  further  shown  to  qualitatively 
predict  the  distribution  of  planar  heat  generation  in  the  laminates.  Experimental  data  from 
laminate  surface  temperature  measurements  using  liquid-crystal  sheets  compared  well 
qualitatively  with  the  theory.  As  a  result  of  this  study,  a  fundamental  understanding  of  the 
controlling  mechanisms  of  thermal  generation  in  continuous-carbon-fiber  systems  under  the 
influence  of  an  alternating  magnetic  field  is  established. 

In  order  to  correlate  this  planar  electric  potential  distribution  to  thermal  generation,  it  is 
necessary  to  model  the  mechanisms  of  field  distribution  through  the  thickness.  The  information 
presented  in  this  work  can  be  refined  by  taking  into  account  the  through-thickness  response  (i.e., 
How  do  the  nodal  voltages  between  fibers  in  adjacent  orthogonal  plies,  obtained  from  the  planar 
grid  model,  interact  with  each  other  to  establish  linear  electric  fields  in  the  interfiber  polymeric 
regions?).  This  requirement  is  accomplished  through  the  fiber  layer  submodel  to  be  presented  in 
a  separate  communication. 

The  planar  grid  model  and  the  supporting  experimental  evidence  address  a  new  complexity 
to  the  issue  of  joining  and  field  repair  of  thermoplastic-based  composites  by  magnetic  induction 
heating.  The  possibility  of  extreme  gradients  of  heat  generation  in  the  plane  of  these  materials 
demands  further  research  in  this  area. 
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